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We explore the possibility that, at zero temperature, vortices can be created spontaneously in 
a condensate of cold Fermi atoms, whose scattering is controlled by a narrow Feshbach resonance, 
by rapid magnetic tuning from the BEC to BCS regime. This could be achievable with current 
experimental techniques. 
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Causality imposes strong bounds on a system whose 
environment is changing rapidly. If external influences, 
such as quenches, try to make correlation lengths change 
faster than the relevant causal speed (e.g. the speed of 
sound) then the system will freeze, to unfreeze later when 
the effect has diminished. As proposed by Kibble and 
Zurek 0-0], this then leads to frustration. In particular, 
vortices will arise spontaneously if they can be trapped 
by domain boundaries. 

Spontaneous vorticity has already been observed in 
cold bose gases, on condensing them by rapid cooling Q. 
In this letter we suggest, on similar causal grounds, that 
vortices can be created spontaneously in a condensate of 
cold fermi atoms by rapidly tuning the binding energy 
of a dominant Feshbach resonance [j| with an external 
magnetic field. 

The idea is very simple. Weak fcrmionic pairing gives 
a BCS theory of Cooper pairs, whereas strong fcrmionic 
pairing gives a BEC theory of diatomic molecules. The 
crossover is characterised by the divergence of the s-wave 
scattering length as @ • On driving the condensate from 
the deep BEC regime to the deep BCS regime by ramping 
an external magnetic field H, the speed of sound v s in- 
creases from essentially zero to 0(vf), the Fermi velocity. 
However, the correlation length £ decreases as the veloc- 
ity increases, from a high value when the initial speed 
of sound is sufficiently small. Thus, if the crossover is 
fast enough the condensate has to be frozen initially to 
prevent the correlation length collapsing acausally fast. 

An estimate of the time t at which the system unfreezes 
is [l|, Q when it can change no faster i.e. 



\i(t)\^ Vs (t). 
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When this happens, with finite correlation length, vor- 
tices will appear to accommodate the frustration of the 
field. In the KZ scenario it is suggested that vortex sepa- 
ration at their time of spontaneous production is £ ~ £©• 
If £ =kp 1 , the inverse Fermi momentum which sets the 
atomic separation scale, and tq = Ti/tF, the inverse Fermi 
energy (in units of h) , we shall show that 



provided tq 3> tq. The timescale tq is the quench time 
for the change in the inverse scattering length induced 
by the changing magnetic field, and is proportional to 
the quench time tjj for the field change. Experimentally, 
with current techniques, tq can be made comparable to 
To itself, suggesting that spontaneous vortex production 
should be observable in realistic systems. In this regard 
there are many similarities with the analysis of [7[ for 
thermal quenching of condensates (although in that case 
causality is determined by second sound) . We stress that 
the mechanism that we are invoking here differs from 
that of defect formation in quantum phase transitions in 
T = condensates . 

For the sake of analytic simplicity, we restrict ourselves 
to narrow Feshbach resonances. Our starting point is the 
exemplary 'two-channel' microscopic action (in units in 
which h = 1) 



S = J dtd s x^^( 
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- g [p(x) ^(x) ^(x) + <Kx)i%(x) 1>l(x)] j (3) 

for cold (T = 0) fermi fields tp a with spin label a = (j, 1), 
which possess a narrow bound-state (Feshbach) reso- 
nance with tunable binding energy v, represented by a 
diatomic field 4> with mass M = 2m. This model has been 
discussed on great detail by Gurarie and Radzihovsky @ 
and we borrow several results from their paper. 

S is quadratic in the fermi fields. Integrating them out 
[Tol | enables us to write S in the non-local form 



Snl = —iTrlnQ 

dt d 3 xcj)*(x) 



d t + 



2M 



2fi — i/ 



<f>(xU) 



6 ~ £o(iq/t ) 
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in which Q 1 is the inverse Nambu Green function, 

c-\ _ f idt-e -g<t>{x) \ , , 

y ~ 1 -9<t>*{x) idt + e J [b) 
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where —g4>(x) — g\(j)(x)\ e lB ^ represents the condensate 
(and e = -V 2 /2to- /i). 

The s-wave scattering length as is determined from 
the binding energy as 



2/i — v 



9 2 m 



(6) 



To see the effect of applying an external magnetic field 
H, we adopt the parametrisation 



as = atg 1 



H-Hn 



whence 



2p 



g 2 m 



<i-Ka bg H b 



(H-H ), 



(7) 
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where a bg is the background (off-resonance) scattering 
length and H u the so-called 'resonance width' @. H 
is the field required to achieve infinite scattering length 
(the unitary limit). As H increases through Hq we pass 
from the BEC to the BCS regimes (with a$ going from 
+vc to -ve). 

In this paper we restrict ourselves to the mean-field 
approximation, the general solution to 5Snl = 0, valid 
if <f> is a sufficiently narrow resonance [t| [ll| . The action 
possesses a U(l) invariance under 8 — > 9+ const., which is 
spontaneously broken. 5 Snl = permits the spacetime 
constant gap solution \(f>(x)\ = |</>o| ^ 0. We perturb 
in the derivatives of 9 and the small fluctuation in the 
condensate density 5\(f>\ — \(f>\ — \4>o\ and its derivatives. 
9(x ) is not small. Using the results of our earlier papers 
[HI IHI; f° r which (|3|) is a limiting case, we can extract 
from Sn l a local effective Lagrangian density L e f / , 

Leff = -\ Po G{9,e) + ^G 2 {9,e) 

-aeG(9, e) + \ijX 2 (e, 9) - -M 2 e 2 , (9) 

valid for long wavelength, low-frequency phenomena. 

L e ff is given in terms of the Galilean scalar combi- 
nations G(9,e) = 9 + (V6») 2 /4m + (Ve) 2 /4m, X(e,6) = 
e + V#.Ve/2m, where the dimensionless e ex S\cj)\ is itself 
a scalar. We postpone a discussion of the coefficients of 
© , except that po = p$ + Pq is the total (fixed) fermion 
number density where p$ is the explicit fermion density, 
and Pq — 2|0o| 2 is due to molecules (two fermions per 
molecule) . For the evolving system the molecular density 
is p B = 2\(j)\ 2 = +4<5|0||0 O |, showing that e cx 5p B , 
the molecular density fluctuation. Although the details 
are immaterial, we have scaled e so that it has the same 
coefficients as 9 in its spatial derivatives. 

For small fluctuations, the linear approximation to the 
Euler-Lagrange equations for 9 and e, sufficient to deter- 
mine the speed of sound, is 
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(10) 



On diagonalising, we see that for long wavelengths the 
gaplcss phonon has dispersion relation u> 2 — w 2 k 2 , 



Po/2m 



N +4a 2 /M 2 



(11) 



independent of ip With 77-dependent coefficients, in the 
deep BCS regime (small a) v s —> vbcs = vf/V^ and in 
the deep BEC regime (large a) v s —> 0. 

Eqs. pT])) represent a two-component system of 
molecules and atom pairs. Consider a spatially homo- 
geneous condensate. Then these linearised EL equations, 
which ignore damping, display the oscillatory solution 



6p B = Spfi cos fit, 



(12) 



describing the repeated dissociation of molecules into 
atom pairs and their reconversion into molecules. The 
frequency Q of density fluctuations can be shown to 
increase monotonically from the exponentially damped 
g\(po\ = 0(/j,exp(—TT/2kF\as\) in the BCS regime to 



ji as it crosses into the BEC regime 



[JJ, [JJ]. In the deep BEC regime Q, 2p. 

A two-component density in which fermions oscillate 
from one to the other is not the same as a two-fluid pic- 
ture. In general, density fluctuations act as sources and 
sinks in the continuity equation and the condensate be- 
haves as a fluid only when these are ignorable i.e. we 
can neglect the spatial and temporal variation of e, in 
comparison to e itself. As a result, e w — 2aG(9) / M 2 , a 
slave to the phase. We shall see that, for our quenches, 
this is the case, whence the Euler-Lagrange equation for 
9 is, indeed, the continuity equation of a single fluid, 



0_ 

dt 



P + v • (H = o, 



(13) 



with p = p a +2ae~N G{9) and v = V9/2m. To complete 
the fluid picture we observe that this definition of p is no 
more than the Bernoulli equation 



1 , 
Sh + —m v' 



= 0, 



(14) 



where the enthalpy 5h = mv 2 5p/p. 

Since the hydrodynamic equations can be derived from 
a Gross-Pitaevskii (GP) equation, it is more transpar- 
ent to reconstitute this equation, with its natural vor- 
tex solutions. Consider the Lagrangian describing the 
wave-function ip of a particle of mass 2to, interacting 
non-lincarly with itself, 



lwo = mr^ - -^vr ■ V0 s -(\ip\ 2 - po) 2 (is) 

4m po 



where we have restored factors of Ti. The Gross-Pitaevskii 
equation following from ((T5|) is 



ihip ■ 



2m 



V 2 ip + 2mv 2 s il) 



2mv\ 
Po 



V#| 2 = 0. (16) 
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FIG. 1: A condensate with g = 1.5, quenched by a mag- 
netic field with quench rate tq/tq = 16. For the solid line 
the ordinate is v^/vp and for the dashed line the ordinate is 
y/h/2M(d/dt)v%. In driving the system from the BEC regime 
(on the right) towards the BCS regime (on the left) the sys- 
tem unfreezes at 1/kpas ~ 0.75, when the two curves diverge, 
from (|18p . The inset shows pQ /po (dashed curve) and Po /po 
(solid curve). 



with coefficients varying smoothly as we cross the 
unitary limit. If wc set ip = ^J~p exp(i#) and solve (fT6|) 
at the relevant order in derivatives, we recover (|13l) 
and |T3|) on ignoring e density fluctuations (see [15| 
for a comparable analysis). Remarkably, the depen- 
dence of the GP equation on the coefficients of © is 
implicit, through v 2 . The length scale of the system is 
£ = h/2mv s = Ti/Mvs and the time scale is r = Ti/Mv 2 . 

The single-fluid model requires that density fluctu- 
ations can be averaged to zero on the timcscalc r i.e. 
rfi ^ 1 which, for the moment, we assume. We stress 
the importance of the narrowness of the resonance. 
At best, broad resonances require a two-fluid descrip- 
tion [10l | from which we have no simple way to draw 
conclusions about defect formation in a field quench. 

Suppose that H increases uniformly in time with 
H/H\h (i = th^ 1 , where t = is the time at which the 
system is at the unitary limit. We write 1/kpas = ~t/ T Q 
for small t, where 



tq = t h 



Ho~ 



(17) 



The time t of ([T]) at which the system unfreezes satisfies 

h d 



{<{t)Y 



2M dt 



«;(*)■ 



(18) 



Empirically, we find the approximate universal behaviour 

v 2 s (t)^v 2 F (ro/ir Q ) (19) 

for a very wide range of parameters. This universality is 
not surprising since, to a fair approximation, v 2 /vp can 
be approximated as 



v l ~ ( w f/ 6 )[1 + tanh(co - b /k F a s )], 



(20) 



(of which we see only the tail in Fig.l) where &o(ff) is ap- 
proximately constant for a wide range of g, which gives 
v s/ v f T o/ T Q automatically. We only need the coeffi- 
cients of L e f f of © to improve upon (|2"0"|) and determine 
the coefficient of proportionality. These can be obtained 
from and [l| ( in the limit of no contact interac- 
tions). The correlation length at the time of unfreezing 
then satisfies the scaling law (T2J) 



^(tq/to) 1 / 2 



(21) 



provided tq ^> tq. 

Eq. (|16p permits vortex solutions within which £ de- 
termines vortex width (in fermion number density). We 
can ignore the order parameter fluctuations 6\ip\, since 
they are shorter-ranged and do not affect how vortices 
pack. The KZ picture then gives an estimated vortex 
separation at time of production as (|2"Tj) . as claimed in 
the introduction. 

The quench parameters are related to the width of the 
resonance r by [gj T m imfi 2 B al g H 2 /h 2 , where [Ib is 
the Bohr magneton. In practice, it is more convenient 
to work with the dimensionless width 70 ~ y/To/ep, 
whereby 



TQ 
TQ 



TQi-F 



7o- 



(22) 



To be concrete, consider the resonance in 6 Li at Hq = 
543.25G, discussed in some detail in |16| . As our bench- 
mark we take the achievable number density po ~ 3 x 



10 1 



whence cf ~ 7 x 10 eV and 70 



0.2. 



In terms of the dimensionless coupling g, where g 2 = 
(64e 2 ? /3/c|.)g 2 , 6 Li at the density above corresponds to 
g 2 < 1. For a condensate of density p it follows that 

f * (Po/p) 1/3 and 



19. 

TO 



H\po 



(23) 



where H is measured in units of Gauss (ms) . Experi- 
mentally, it is possible to achieve quench rates as fast as 
H w O.lG/ms [Hj]. 

The condition rfl ^> 1 throughout the quench now 
becomes 



n 



rfi = — 



ep 4vJ 



Q Tq 

— ^ > 1. 

e-F to 



(24) 



With il sa 2\fj,\, large in the BEC regime, the inequality is 
guaranteed there for all tq /tq 1. However, in the deep 
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BCS regime Q is exponentially damped, not sustaining 
the single fluid picture. Thus the single fluid approxi- 
mation, on which our causal analysis depends, requires 
that the quench does not trespass beyond the unitarity 
limit. The question is then whether such a quench can be 
implemented sufficiently fast for the condensate to freeze. 

We believe it possible for realistic quenches. It is sim- 
plest to adopt a density p < po, permitting smaller H to 
maintain tq > tq. As an example (sec Fig.l), we take 
g = 1.5, corresponding to p w O.lpo, and tq/t = 16. 
Then, the domain structure is formed when v s rj vbcs/L 
with 1/kpas ~ 0.75 and Pq w 0.38po,pQ rj 0.62/?o- tH 
decreases throughout the quench, but its final value of 
16f2/eF ~ 27 remains sufficiently large to justify the ap- 
proximation. Adjacent parameter values are equally suc- 
cessful. 

This suggests that, with £ ss 4£o, spontaneous vortex 
creation should be possible, since the length scale £ c for 
a condensate of N = 10 5 atoms at this density would 
give £ c w 100£o- We should beware that the extent to 
which the KZ bound ^ is saturated depends on the sys- 
tem. To cite extremes for spontaneous vortex formation 
at thermal quenches, it is saturated for vortex produc- 
tion on quenching 3 He — B 17 1, but underestimates vor- 
tex separation strongly for high-T c superconductors [l8[ . 
The spontaneous creation of vortices in thermal quenches 
on low-T c superconductors [l9| and fluxons in Josephson 
junctions |20j give results in between. Finally, we see that 
the approximation from Eq.([7|) to Eq.® is justified, with 
a fractional error less than kpdb g <C 2.10~ 2 . 

Our conclusions that spontaneous vortex creation 
could be achievable with current experimental techniques 
needs further qualification. We stress that, since the vor- 
tices form early in the ramp we do not have to continue 
it far into the BCS regime, nor begin it from the deep 
BEC regime, where our our idealised narrow-resonance 
approximation fails. This is rather like the situation in 



thermal quenches in which defects form so close to the 
critical temperature that there is no need to cool much 
below it. This has the further advantage in that, al- 
though our idealised calculations were for temperature T 
= 0, in reality temperature is finite. By stopping soon 
enough, we would hope to remain clear of critical thermal 
behaviour. 

Narrow resonances are difficult to work with because 
of the required field stability, but we expect them to give 
most defects after a ramp. Increasing resonance width 
in (f2"2"|) increases tq and hence £ at fixed density. How- 
ever, with £ cx 7 for moderately narrow resonances, 
the effect of broadening the resonance is, initially, weak 
and we can still anticipate observable spontaneous phase 
change for large condensates. [For very broad resonances 
we have no reliable analytic causal constraint for £.] 

As a final caveat we do not have the homogeneous 
condensates assumed above and should take the de- 
tails of their trapping into account. The causal length 
£ cx (poa/To) 1 ^ 2 depends upon density and will vary 
across the trap, but vortices should still form in its centre 
if the condensate is sufficiently large, albeit with profile- 
dependent allometric scaling behaviour. In this regard 
there are many similarities with the analysis of [7j for 
thermal condensates and we would have to modify our 
analysis appropriately. This letter is rather aiming for a 
proof of principle, that causality could lead to observable 
changes of phase accessible by current experiments. 

RR would like to thank the National Dong Hwa Uni- 
versity, Hua-Lien, for support and hospitality, where 
much of this work was performed and Profs. Randy Hulct 
and Matt Davis for helpful comments. The work of DSL 
and CYL was supported in part by the National Science 
Council and the National Center for Theoretical Sciences, 
Taiwan. 



[1] T.W.B. Kibble, Physics Reports 67, 183 (1980). 

[2] W.H. Zurek, Nature 317, 505 (1985), Acta Physica 

Polomca B 24, 1301 (1993). 
[3] W.H. Zurek, Physics Reports 276, 4 (1996). 
[4] C.N. Weiler et al., Nature 455, 948 (2008) 
[5] M. Greiner, C. A. Regal, and D. S. Jin, Nature 426, 

537 (2003); S. Jochim et al. Science 302 (2003); M. W. 

Zwierlein et al., Phys. Rev. Lett. 91 250401 (2003). 
[6] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 

92, 040403 (2004); M. W. Zwierlein et al, Phys. Rev. 

Lett. 92, 120403 (2004); C. Chin et al, Science 305, 

1128 (2004); Y. Shin et al., Nature (London) 451, 689 

(2008). 

[7] W. H. Zurek, Phys. Rev. Lett. 102, 105702 (2009). 
[8] B. Damski and W.H. Zurek, Phys. Rev. Lett. 99, 130402 
(2007). 

[9] V. Gurarie and L. Radzihovsky, Annals Phys. 322, 2 
(2007). 

[10] D-S. Lee, C-Y. Lin and R. J. Rivers, Phys. Rev. Lett. 98, 
020603 (2007) and references therein. 



[11] S. Diehl and C. Wetterich, Phys. Rev. A 73, 033615 
(2006). 

[12] C-Y. Lin, D-S. Lee and R. J. Rivers, Phys. Rev. A 80, 

043621 (2009). 
[13] V. Gurarie, Phys. Rev.Lett. 103, 075301 (2009). 
[14] E. Timmermans, K. Furuya, P. W. Milonni, A. K. Ker- 

man, Phys. Lett. A 285 228 (2001). 
[15] I. J. R. Aitchison, P. Ao, D. J. Thouless and X.-M. Zhu, 

Phys. Rev. B 51, 6531 (1995). 
[16] K. E. Strecker, G. B. Partridge and R. G. Hulet, Phys. 

Rev. Lett. 91 080406 (2003) 
[17] V.M.H. Ruutu et al, Nature 382, 334 (1996). 
[18] A. Maniv, E. Polturak, G. Koren, Phys. Rev. Lett. 91, 

197001 (2003). 

[19] R. Monaco, J. Mygind, V.P. Koshelets and R.J. Rivers, 

Phys. Rev. B80 (2009) 180501 
[20] R. Monaco, J. Mygind, M. Aaroe, V.P. Koshelets and 

R.J. Rivers, Phys. Rev. Lett. 96, 180604 (2006) 



